Exact solutions of the two-mode model of multicomponent Bose-Einstein condensates 
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We find the explicit solution of a two-mode model used to explain vortex dynamics in multicompo- 
nent Bose-Einstein condensates. We prove that all the solutions are constants or periodic functions 
and give explicit formulae for the time evolution of the populations of the two atomic species present 
in the condensate. 
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I. INTRODUCTION 

One of the most remarkable achievements of contemporary 
physics has been the realization of Bose-Einstein condensa- 
tion with ultracold atomic gases In those experiments 
bosonic neutral gases where cooled down below the critical 
temperature and a collective coherent behavior was observed 
in the gas cloud. 

One later interesting results in this field were the achieve- 
ment of simultaneous condensation of several atomic species 
using sympathetic cooling Q and the generation of multi- 
ple condensates using r.f. transitions between different hy- 
perfine levels of Rb*^ ||. These two species are usually de- 
noted by |1) and |2). Despite its complexity, this system may 
be described in the mean field limit using two macroscopic 
wavefunctions describing each of the atomic species. After 
these achievements many theoretical and experimental works 
on multicomponent systems followed 

The first generation of vortices in Bose-Einstein conden- 
sates reported in Ref . jsj , has been one of the most remarkable 
goals which has been achieved using multicomponent systems. 
From Ref. we know that while one may build two possible 
configurations with a unit charge vortex only one of them is 
stable. The stable configuration corresponds to the vortex 
placed on the |1) state, namely the one with the largest scat- 
tering length. The configuration with the vortex placed in the 
|2) state, on the other hand, leads to some kind of instability. 

In recent works numerical simulations have been used 
to show that the dynamics of the unstable configuration can 
be understood within the framework of mean field theories 
for the double-condensate system. In a later work Q, it was 
shown that a simple two-mode model is enough to describe 
accurately the dynamics of the multiple-condensate system in 
two dimensional setups (e.g. oblate condensates where the z- 
coordinate may be neglected). The model was able to describe 
the most relevant features found in experiment and also in the 
mean field theories of the system j^Q. 

In this paper we concentrate on the analysis of the simple 
model reported in Ref. proving that the system is inte- 
grable and all the solutions are periodic or constants (equilib- 
ria). This proves that the tendency of the system to exhibit 
periodic exchange of vortices between both atomic species is 
fundamental and not the consequence of a fortunate choice of 
parameters. We also study the period of the vortex trasfer 
mechanism as a function of the physical parameters of the 
problem, in particular we find a linear dependence on the 
total number of particles A'^ of the system. 



II. TWO MODE MODEL FOR THE MULTIPLE 
SPECIES CONDENSATE. 

A. Gross-Pitaevskii equations for the two-species 
system. 

In this work we will use the zero temperature approxi- 
mation, in which collisions between the condensed and non 
condensed atomic clouds are neglected. In the two species 
case this leads to a pair of coupled Gross-Pitaevskii equations 
(GPE) for the condensate wavefunctions of each species 
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where Uij — A-Kh^aij /m are constants controlling the nonlin- 
ear behavior, which are proportional to the the s-wave scat- 
tering lengths of 1-1, (an), 2-2 (022), and 1-2 (012) binary 
collisions. 

To simplify the formalism we assume that both trapping 
potentials are concentric and spherically symmetric, Vi(r) = 
y2(r) — ^mLpr"^ , just like in the experiment. The consider- 
ation of the differences between Vi and V2 does not add new 
physics to the model as discussed in Ref. 

Next we change to a new s et of u nits based on the trap 
characteristic length, ao = ^Jh/muj, and period, t = 1/uj 
defined as a:: — > x/ao,t — > t/r, Uij a,ij I ao and j (x) — 

Nj'ijjj{'x.). Equations (^) conserve the number of particles on 
each hyperfine level and so we may choose 
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This choice implies that the particle number of each species 
appears in the nonlinear coefficients Uij. 

With the previous rescaling, the GPE for the multicompo- 
nent system read 



— - + -r^ -I- + U\l\'4}2^ 
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The coefficients of the matrix of nonlinear coefficients 
satisfy the relations uwjuxi = aiiA''i/ai2 A2, U21/M22 = 
a\2Ni/a22N2, which means that except for the particular case 
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in which Ni — N2 — N, this matrix is nonsymmetric. In 
terms of the population imbalance l3 — N2 /Ni , and for a fixed 
total number of particles the values are un = gaii/{l + 13), 
U12 = gai2l3/{l + l3), U21 = gai2/{l + f3), U22 = 3a22/3/(l + 
where g = inN/ao. 



B. Derivation of a two-mode model 

For completeness we will derive here the two-mode model 
following a different formalism than that of Ref . . The idea 
is to use a limited expansion 



i>nit,r) = auk{t)<f)k{r) , 



(4) 



where we take 4'k (r) to be the lowest energy solutions of the 
linear problem 
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A + l/(r) 



(/>fc(r) = Ek(f>k{r)- 



(5) 



Although we have choosen a harmonic oscillator basis, any 
suitable truncated basis including a mode resembling the 
ground state and another one corresponding to a state with a 
vortex could be used in this expansion. Despite its simplicity 
and crudeness (only two modes of the full expansion are kept), 
this approximation has been shown to retain many qualita- 
tive features of the dynamics of the system [^. Inserting Eq. 
(^) with (f) given by (g) into Eq. we get 

2 

^ [-ih^k{t) + Eka^k{t) + Un{t, r)a„fc(f)] (j)k{r) = 0, (6) 
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where 



— 2Ci2iii2 P11P12P21P22 sin$, 
+2C12M12 P11P12P21P22 sin$, 
+2C12M21 P11P12P21P22 sin$, 
-2Ci2Ji2i P11P12P21P22 sin$, 



flnfc = Pnfc exp [lH„k ) 



and 



being 
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^nkPnk ? 



711 = Cll (Uii - M21) + C12 {U21 - 2Uii) . 

712 = C22 ("21 - Mil) + C12 (2uii - U21) , 
722 = C'22 (m22 — U12) + C12 {ui2 — 2U22) ■ 
721 = Cll {ui2 - U22) + C12 {2U22 — U12) . 
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These equations must be complemented with the appropriate 
initial conditions for the densities pki and phase $. 



Pnk{to) 
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n,fc = 1,2 
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The form of the two-mode model [Eqs. (p^j^jis])] presented 
here is slightly different from that of Ref. fsTbut is more 
suitable for our purposes in this paper. 



where Un{t,r) — J2m 

\i>mit,r}\\ After some algebra 

we find 

dt \a„kf = — 2u„fiCi2 I aiiai2a2ia22 | sin<l>„fc, (7) 

where Cjk = {\<P]\''' \<Pkf) and $ = $11 = -$12 = -'l>2i = 
$22 = arg(aiia22/ai2C[2i). The equations for the phases are 



dt arg a„fc = -Ek - Re (^fe [/„ (f)k') ■ 



(8) 



where {A) — J d^fA. Calculating (^(pt Un (pk'} for the two 
cases k' — k and k' — k, where k is the complementary value 
of the index, k, k' = 1,2 
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((/ifc Un (j)k) = ^ Unm {Cfcfc \amkf + Cl2 Idmfcl^} , (9) 
m — 1 

2 

((j!>fc Un (l>k) = Cl2 ^ Unm Omfc flmfe, (10) 
m — 1 

we get the following set of equations for our problem 



III. SOLUTION OF THE EQUATIONS. 

Let us first note that the two-mode model equations have 
the following constants of motion: 



M = pn + p?2 = const = p?2 + Pi2, 
^2 = P21 + P22 = const = P21 + P22, 

L = M21P12 + W12P22 = W2lP?2 + W12P22- 
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It follows from (0-(|l|) and ^ that 

ftp = 27,Co2 P11P12P24P22 sin (26) 

where 

7* = Mi2 (712 — 711) + "21 (711 - 722) . (27) 
The following relations are then evident 

(28) 
(29) 
(30) 
(31) 
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with f — 'Y^^f.^nkpiik- These equations imply that all the 
relevant densities may be obtained from a single quantity, 
T{t). Our goal now will be to find an equation ruling its 
dynamics. Let us notice that Eq. (If fl) can be rewritten as 



dt piipi2p2ip22 cos $ = — r dt r 



(32) 



and solved: 



p\\pl2P2ip22 COS $ — p\\pi2p21p22 COS $ 



4Ci27* 



(33) 

Using the last equation and (p8|)-(|3l|), Eq. (^ ) for Vt can be 
presented as 



r? =p4(r) 



(34) 



where 



Vi{T) = ~- [r^ + 47,Ci2pii/5i2/52i/522 Cosl']% 



4C 



2 M?2ll21 
12 ^,2 



(r - Tn) (r - ri2) (r - r2i) (r - r22) (35) 



and the constants Tnk, n,k — 1,2 are defined as Fn = 

f + 7,/Ml2Pn, ri2 = f - 7,/ui2 /5i2, r2i = f - 7,/m21 pil, 

r22 = r + 7,/u2i /322- Thus, the solution of our equations can 
be obtained from a single equation for F, which, using the 
previous equations, allows to obtain the expressions for the 
densities pjk(i). In particular Eq. (M) leads to 
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Let us rewrite the polinomial (^5|) in terms of its roots as 
follows 



P4(F) (F-Pi)(F-P2)(F 



0(r- 



(37) 



where = 4C?2 v.hu2i/'y* - 1/4 and Pi, P2, Qi, (52 are the 
roots of P4 . There are at least two real roots Pi , P2 of P4 
satisfying that F is always located between them, i.e. Pi < 
f < P2. 

To prove the last affirmation let us first compute 
P4(f) = 4C?27; sin^ $ > 0. 
Next we evaluate 

P4(Fii) = -i [F?i -f' +47.Cl2Pll/5l2/521/522 COS<l]^ < 0, 
P4(Fi2) = - i [FL -f^ + 47.Cl2Pll/5l2P2lP22COS<l]^ < 0. 

This means that there is at least one real root located in each 
of the intervals Pi e [Fii,f] and P2 £ [f,Fi2]. 

As a corollary we get that p4(F) > for Pi < F < P2. 
Concerning the other two roots Qi,2 {Q2 > Qi if real), noth- 
ing may be said in general. However, for the case P* > 
it is clear that p4(±oo) > and then Qi,2 are real num- 
bers. When Pt < it happens that Qi, Q2 could be (and, in 
fact, they are) complex conjugate roots for certain parameter 



ranges. In any case, the boundedness of F G [Pi,P2] ensures 
the periodicity of the solutions. 

Let us now proceed to find an explicit form for the solutions. 
First, to present polynomial P4(F) in the canonical form we 
make the transform: 



m 



aY + 13 
Y + 5 ' 



where 



a=i('5 + l)P2-i(5-l)-Pi, 
P=^^{& + l)P2 + ^{5-l)P^, 



5 = 



v'IPi - Qi| |Pi - O2I + ^\P2 - Oil IP2 



^\Pi - Qi\ \Pi - O2I - ^\P2 - Qil IP2 - ' 

Equivalently we may also write 

A Y 

T{t) = Po + ^ 



5 5 + Y 

with the quantities A and Po being given by 
A = a5-/3=i((5"-l) (P2-P1), 



Po 
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[{5 -f f ) P2 + (5 - I) Pi] 



(38) 
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It is clear that Pi < Po < P2 and thus P(Po) > 0. In terms 
of Y the polynomial P4(F) can be written as 



P(F) = <5^P(Po; 



{l-Y^) (l-fc^ya) 



{5 + Yf 

Thus, the function Y{t) satisfies the equation 



Yt=Q, ^(l-y2)(i_fc2y2) 



(45) 



(46) 



where 



y/\Pl - Qi| IP2 - Q2\ + y/\Pi - Q2I IP2 - Ql 
Equation (|4^) can be solved as follows: 

Y{t) = sn (h(t-i) , (48) 

where 



^|Pl - Oil IP2 -Q2I - ^|Pl-Q2||P2-Ql 



\/\Pl-Ql\\P2-Q2\ + y/\Pl-Q2\\P2-Ql\ 



(49) 



and the constant t may be determined from the initial condi- 
tions by solving the equations 



(^( 



to-t) ,k] = 



(50) 



This allows us to get a explicit solution for our problem. To 
do so let us first use Eq. mm to find 



F(t) 



(^( 



t-i) ,k] +p 



t-t) ,k] +5 



(51) 



(47) 
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Thus, inserting this explicit expression for r{t) and Eqs. ([2^)- 
one may get explicit expressions for the populations 
Pjk{t),j,k = 1,2. Typical forms of the solutions for different 
parameter values are shown in Fig. |l|. 
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FIG. 1. Evolution of the densites pii{t) (solid) and p\2{t) 
(dashed) for two set of combinations of parameters and 
initial conditions (a) g = 100, an — 1, ai2 — 0.97, 
022 = 0.94, p = = 1, pii(to) = 0.78, (Piiito) = 0, 

0i2(to) = 0.5,p2i(to) = 0.5, <^2i(fo) = 1.5,(?!)22(fo) = 0.7. (b) 
g = 100, an = l,ai2 = 0.1, 022 = 0.02, /3 = N1/N2 = 1/4, 
pii(to) = 0.9, (^ii(to) = 0, (;!>i2(to) = 0, p2i(to) = 0.9, 
4>2i{'to) ~ 1-5, (/f>22(to) = 0. We fix t = and get the cor- 
responding (irrelevant) value for to from Eq. (po[). 

The period of oscillations of the populations Pjkit) is that 
of the sn function, i.e. 



T = 



4^K{k) _ AAK{k) _ 8K{k) 



n 



f VIA -Q2I |P: 

p/2 

Jo 



(52) 



1 — fc2 sin S is the 



V|Pi -Qil IP2 

where the quantity K{k) 
complete elliptic integral of first kind llOl, 

Concerning the period, the dependence on the physical pa- 
rameter A'^ is easy to obtain. Let us define the new con- 
stants (not dependent on TV") yij — yij/N, Uij = Uij/N, 



/N, 7* = j*/N^. Using these definitions we 
find that the roots Pi, P2,Qi,Q2 scale linearly on A'^. Thus 
k, S and do not depend on and the period scales as 
T{k,N) =N-'^T{k,l). 

Concerning the instability mechanisms described in Refs. 
they are now completely understood within the frame- 
work of Eq. (|l|). For instance in Fig. ^ we plot typical 
solutions in the regime of vortex trasfer for parameter values 
typical of the experiments |^ . Both the stability of the vortex 
when it is placed in |1) and its instability when put on |2) are 
well described by the exact solutions of the two-mode model. 
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FIC. 2. Evolution of the densites pi2{t) and p22(t) (con- 
taining the vorticity of each component) for two set of com- 
binations of initial conditions corresponding respectively to 
(a) a stable vortex in |1 > and (b) an unstable vortex in 
2 >. Parameter values are g — 100, /3 — N1/N2 — 1, 
ail ~ 1, ai2 = 0.97, a22 ~ 0.94. Initial conditions are (a) 
pii(to) = 0.001, (fiiito) = 0, 0i2(to) = 0.1, p2i(to) = 0.999, 
<l>2i{to) = 0, <?!.22(0) = 0.2.(b) pii(fo) = 0.9999, <?!>ii(to) = 0, 
Mto) = 0.1, p2i(to) = 0.001, <l>2i{to) = 0, <^22(to) = 0.2. We 
fix f = and get the corresponding (irrelevant) value for to 
from Eq. M). 
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IV. CONCLUSIONS 



In this paper we have integrated the two-mode model de- 
veloped for the explanation of the vortex transfer mechanisms 
described in Refs. The most remarkable result is that 

all the solutions are periodic functions and thus the transfer 
mechanisms of at least part of the vorticity are natural within 
the range of validity of the model, which was previously found 
to be at least that of pancake traps and some regimes of fully 
three-dimensional traps. It has been shown that the frequency 
of the oscillations depends linearly on the total number of 
particles of the condensate while the dependence on the other 
parameters (relative populations and scattering lengths) is 
nontrivial and given by our explicit formulae. 

We hope that the technique described here will be useful 
to analyze the similar problem which arises when a Joseph- 
son coupling between both species is incorporated into the 
experimental setup by means of a off-resonant laser |n| . 
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